A climatic thermostat making Earth habitable 
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atmospheres is determined by the radiative balance between the 
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q . non-reflected incoming solar radiation and the outgoing long-wave 
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^ ; black-body radiation from the atmosphere. The surface tempera- 

> : 

•'-j ■ ture is higher than the black-body temperature due to the green- 

house warming. Balancing the ice-albedo cooling and the greenhouse 
warming gives rise to two stable climate states. A cold climate state 
with a completely ice-covered planet, called Snowball Earth, and a 
warm state similar to our present climate where greenhouse warming 
prevents the total glacition. The warm state has dominated Earth in 
most of its geological history despite a 30 % fainter young Sun. The 
warming could have been controlled by a greenhouse thermostat op- 
erating by temperature control of the weathering process depleting 
the atmosphere from CO2. This temperature control has permitted 



life to evolve as early as the end of the heavy bombartment 4 billion 
years ago. 



Introduction 

Primitive life has existed on Earth since early in its geological history. Stro- 
matolites, which are thought to be fosils of photosythesizing cyanobacteria are 
found as old as 3.5 Ga (giga-anni = billion years) QSchopf, 1993| ). Isotope frac- 
tionation in carbon found in 3.8 Ga old rocks form Isua, Greenland indicates 



a biological origin ( Rosing, 1999 ). This is perhaps within a few hundred mil- 
lion years the earlies time of stable planetary climate possible for life. In the 
prior heavy-bombardment period impacts releasing energies enough to evapo- 
rate the entire ocean would probably have sterilized the Earth if life existed at 
that time (|Kasting fc Catling, 2003). 



At some time around the Archaean-Proterozoic border, 2.5 Ga BP, the 
atmospheric content of oxygen began rising to its present level, making way 
for dominance of aerobic life forms. For a very long period after then there is 
no evidence of biological evolution. The eukaryotic (cells with a nucleus) and 
multicellular life seems to have evolved as late as 0.6 Ga BP at the Cambrian 
explosion, or perhaps in the late Precambrium where the soft body Ediacara 
fauna originates. This could indicate that the oxygen level rose slowly, reaching 
levels needed for effective oxygen transports and metabolisms in multicellular 
organisms at the Cambrian explosion. 

The habitable zone (HZ) around a sun-like star is defined by the require- 
ment of liquid water at the surface of planets or moons in the zone. This 



zone is primarily determined by the radiative output from the star. The de- 
pendence of the radiative flux, which goes as inverse distance squared makes 
the HZ relatively narrow. At the present time our solar system contains only 
Earth within the HZ flKasting fc Catling, 2003). 



The long existence of life indicates that the surface temperature on Earth 
has indeed kept within the narrow range permitting liquid water for sur- 
prisingly long, despite a 30 % lower solar flux at the beginning of Earth's 
history. This apparent paradox is dobted "The Faint Young Sun Problem" 
(|Kasting, 19880. The solution to the problem could be a selfregulatory mech- 



anism, either through a biotic feedback ( Lovelock fc Margulis, 1982 1 or a geo 



chemical regulation ( Owen et al, 1979] Walker et al, 1981 ). The former can 



only be at play after the evolution and expansion of life to the planetary scale, 
and can therefore not explain the stable climate conditions necessary for the 
initiation of life. 

A geochemical regulation could, besides extending the time for habitability 
of a planet, potentially widen the HZ. There is now mouting evidence that Mars 
had liquid water on the surface in its early history and was thus within the 
HZ then despite the lower solar luminosity. 



The radiative energy balance 

The surface temperature of the planet is determined from the balance between 
incoming solar radiation Ri and outgoing black-body radiation R . This energy 
balance depends, among other factors, on the planetary albedo. The albedo of 
an object is the fraction of the sunlight hitting the object which is reflected. 



The planetary albedo is not a constant factor; it depends, through the amount 
of clouds and ice, on the state of the climate itself. The feedback of clouds on 
temperature is very complicated. It depends on the height in the atmosphere 
where the clouds are and the state of the atmosphere surrounding the clouds. 
The clouds cool by reflecting the incoming radiation and they heat by trapping 
the outgoing radiation. Ice and snow on the surface unambiguously cool by 
reflecting the incoming short-wave radiation, so the amount of ice and snow 
influences the planetary albedo. This effect can be described in a model of 
the climate represented by just one parameter, the mean surface temperature 
T ( Crafoord fc Kallen, 1978| Ghil et al, 19 85). This temperature determines 



the long-wave radiation and the reflection of the short-wave radiation through 
the albedo. The amount of ice and snow is larger when the temperature is 
lower so the lower the temperature the higher albedo. If the temperature is 
below some low temperature 7\ the planet will be completely ice covered and 
a further decrease in temperature cannot increase the albedo above the value 
a± for the fully ice covered planet. If the temperature is above some other 
high temperature T 2 the ice is completely melted and a further increase in 
temperature will not lead to a decrease in albedo below the value a-i for the ice 
free planet. The simplest functional form is a linear dependence of the albedo 
on temperature in between these two temperatures. This is a reasonable choice 
when no other information is available a priory. We then have the relation 

(rr\ t (?2 - T)a>i + (T - Ti)a 2 

J 2 — J-l 

where T( a ,b] is the indicator function for the interval (a, b}. The change of 
the temperature T is determined by the difference R\ — R Q in incoming and 



outgoing radiation according to the equation 

= Ri - R a = [1 - a{T)]S - ag(T)T* (2) 

where c is the heat capacity, a is the Stefan-Boltzmann constant and 5 = 5/4 
is a quarter of the solar constant. (The quarter comes from the ratio of the 
cross-sectional area to the surface area of the sphere). The factor g(T) ex- 
presses the atmospheric greenhouse effect. The black-body temperature 
of the planet is the temperature at the level in the atmosphere from where the 
long-wave radiation is emitted. This level is the height of the optical thickness 
at the long-wave band seen from space. Depending on greenhouse gasses and 
clouds the level of outgoing radiation is approximately 3 kilometers above the 
surface. The difference between the black-body temperature and the surface 
temperature is the greenhouse warming (or cooling). On Earth the atmosphere 
is transparent to the sunlight which thus heats the surface, this, in turn, heats 
the atmosphere from below. The lower atmosphere (the troposhere) thus expe- 
riences a negative lapserate (temperature change with height). The lapserate 
depends in a complicated way on the static stability and atmospheric dynam- 
ics. In the present climatic conditions the lapserate is of the order -lOK/km, 
thus the greenhouse effect on Earth is approximately 3kmx 10K/km=30K. 
Without the greenhouse effect there would be no liquid water at the surface 
of the Earth. The atmospheric greenhouse effect, the change in cloudiness 
and other factors must all be expressed through the "transfer function" g(T), 
where T represents a mean surface temperature, from which the black-body 
temperature is derived; ag(T)T 4 = aT^ ^. The simple model (J2J) is relevant 
because the behavior does not depend critically on the specific choice of the 



parameters QBudyko, 1969[ |Sellers, 1969[ |Crafoord fc Kallen, 1978 1 . 



The behavior of (J2J) is easily understood from a graphic representation. 
Figure ^ shows the incoming and outgoing radiations as a function of temper- 
ature. There are three temperatures T a ,Tb,T c for which the curves cross such 
that the incoming and outgoing radiations are in balance. These points, the 
fixed points, are the stationary solutions to (J2J). Consider the climate to be 
at point T a . If some small perturbation makes the temperature become lower 
than T a we will have Ri > R D implying that cdT/dt > and the temperature 
will rise to T a (see figure 1). If, on the other hand, the perturbation is positive 
and the temperature is a small amount larger than T a we have Ri < R a imply- 
ing that cdT/dt < and the temperature will decrease to T a again. Thus T a is 
a stable fixed point. A similar analysis shows that T& is an unstable fixed point 
and T c is a stable fixed point. If the temperature at some initial time is lower 
than T& it will eventually reach the temperature T a and if it is higher than T& 
it will reach the temperature T c . The present climate is the climate state T c 
where the ice albedo does not play a significant role in cooling the Earth. 

The outgoing radiation depends on the surface temperature T through the 
atmospheric concentration of greenhouse gasses and clouds. Expressing the 
greenhouse effect in terms of the difference between the the surface tempera- 
ture and the black-body temperature T g — (T — the outgoing radiation 
may be written g(T)aT A = a(T - T g - 0.3 * [T — 270]) 4 = R (T, T g ). The 
last factor in the paranthesis is an empirical factor expressing the increase in 
greenhouse effect with temperature due to increase of atmospheric water vapor 
with temperature. 

Consider the climate state represented by T c in the situation T g < T g [Present day] . 



Then for T g not too small, corresponding to the upper dashed blue curve in 
figure HI the equilibrium temperature is lowered a little, as we would expect 
when the greenhouse warming decreases. However, if T g becomes smaller than 
some value Tg the two curves do not cross in more that one point and there 
is no stable fixed point near T c . The climate will then run into the only stable 
fixed point T a which is still present. A saddle- node bifurcation has occurred 
resulting in a large change of climate. If T g grows again the climate state T c 
will not recover until T g exceeds some other value T g 2 ^ > T g [Present day] and 
the system returns through a hysteresis loop. For each value of T g we have 
either one or three fixed points and we can plot the fixed points in a bifurcation 
diagram as functions of T g (figure EJ). The two full curves represent the stable 
fixed points and the middle dashed curve represents the unstable fixed point. 
The unstable and one of the stable points coincide at the bifurcation points 
T« and Tj 2) . 

The stable climate state T a corresponds to a totally ice covered planet. The 
totally ice covered planet has been termed " Snowball Earth" ( Hoffman et al, 1998 ) 



There is geological evidence of such an extreme "deep freeze" climate several 
times in the late Neoproterozoic period around 0.7 Ga BP. This is based on 
findings of glacial deposits like moraine in many places which at those times 
were near the equator. The speculated way out of the deep freeze is the follow- 
ing: The balance in geological timescales between silicate weathering, binding 
atmospheric CO2 into rocks and volcanic out-gassing of CO2 was changed dur- 
ing the deep freeze. Due to the cold conditions the atmosphere dried out and 
weathering effectively stopped. Unchanged volcanic out-gassing resulted in an 
ever increasing amount of CO2 in the atmosphere. At some point, after about 



80 million years, this would result in a greenhouse warming strong enough to 
melt the ice. The warming would then be almost explosive with global mean 
temperature going from some — 40°C to some +50°C within a few years. This 
kind of dramatic climatic changes will strongly stress the planetary biota. 



Stability of the climatic state 

The planetary climate is influenced by internal and external factors which can 
push the climate state away from the equilibrium position. If the perturbations 
are small the equilibrium state will be restored within a typical timescale, 
depending on the size of the perturbation and the strength of the restoring 
force. The fluctuations can be represented as an independent white noise rj (t) 
with intensity a. Linearizing (J2J) around the equilibrium state T c we get: 



c— = -a(T - T c ) + or). (3) 
The parameter —a is the expansion coefficient for the right hand side of (J2J). 



Equation (j3J) is the Ohrnstein-Uhlenbeck process ( Gardiner, 1985 ), where the 
variance is ((T — T c ) 2 ) = a 2 /2a (see figure 3) and the timescale for restoring 
the equilibrium temperature can be defined as the autocorrelation time r = 
oT x . The stability against the random perturbations is thus measured by 
a, the larger a the smaller is the response to the "noise" and the faster the 
perturbation is forgotten. The stability of the climate state represented by the 
temperature T c governs the conditions for biota. However, a long timescale 
stability is not ensured if the parameters governing the value of T c itself is 
changing. The early sun was about 30 % less luminous than today, which 



implies that T c determined by (0) would be lower than today. For a fixed 
value of the greenhouse gas mixing ratio the equilibrium temperature can be 
found from (J2J) as a function of the solar flux S. This is shown schematically 
in figure EJ At times prior to approximately 2 Ga before present the solar 
flux was lower and not permitting liquid water at the surface of the Earth, 
contrary to what is observed, except for the spells of Snowball Earth in the 
late Neoproterozoic period (0.7-0.6 Ga BP). A possible solution to this enigma 
is that the concentration of greenhouse gasses was much higher in the early 
atmosphere. 



The greenhouse thermostat 

The atmospheric concentration of CO2 in geological timescales depends on the 
balance between outgassing through volcanos and burrial through weathering 
of silicate rocks. The rates of change of atmospheric CO2 concentration gov- 
erned by these two factors does not depend on the concentration itself. This 
implies that they will not balance unless they are exactly equal and opposite. 
The outgassing is independent of the surface temperature while the silicate 
weathering rate is stongly depending on temperature. Weathering requires 
precipitation dissolving atmospheric CO2 in the form of carbonic acid. The 
rate of precipitation is stongly temperature dependent. In the present climatic 
conditions weathering takes place in the tropics while it is almost absent in 
the dry and cold polar regions. The simplest way to describe the effect of 



temperature on silicate weathering is by a step function; 



&A = F o -F w 9(T-T ) (4) 

where F Q is the rate of outgassing, F w is the rate of weathering, 9(T — T ) is 
the heaviside step function and T = 285K is the temperature below which the 



is no weathering Q Walker et al, 1981| ). If F w > F Q the CO2 will be depleted 
from the atmosphere. This will, however, cool the planet by diminishing the 
greenhouse effect. When the temperature falls to T weathering stops and the 
atmospheric CO2 rises again by outgassing. In terms of the energy balance 
the greenhouse factor T g should represent the part of the greenhouse which is 
independent of time while the CO2 greenhouse effect must be described as an 
additional factor. The energy balance (J2J) then becomes: 

c— = [1 - a(T)]S - a(T -T g - 0.3[T - 270]) 4 + fO(T - T) (5) 

where f9(To — T) is the additional CO2 greenhouse heating. The behavior 
is now fundamentally different from the behavior of (J2J). For T c < Tq the warm 
climate state is given by 

T< = [{1 - a 2 )S + //a} 1 / 4 + T 9 - 51]/0.8 (6) 

If Tf > T we have 

T> = [{1 - a 2 )S/a} 1 / i + T g - 51]/0.8 (7) 

But T> < Tf =^> T< < T contrary to the assumption. In this case T c is not 
a steady state solution. The surface temperature will increase until it reaches 
To where weathering efficiently depletes the atmosphere from CO2 and the 



temperature drops again. This mechanism is a greenhouse thermostat. The 
weaker heating from the faint young sun is thus automatically compensated by 
a stronger early CO2 greenhouse. The resulting constant surface temperature 
is indicated as the fat blue curve in figure 4. In the future the thermostat 
will not be functional since the atmosphere will be depleted from CO2 so the 
increasing solar flux will result in a steady increase in surface temperature. 



Perspective 

The surface temperature on Earth has been relatively stable permitting liquid 
water necessary for biology as we know it in most of its geological history. 
This is despite a relatively large change in solar flux and atmospheric con- 
tent of greenhouse gasses which suggests a thermostat mechanism controlling 
the surface temperature. A possible thermostat could be operational through 
silicate weathering depleting the atmosphere from the constantly outgassed 
volcanic CO2. The temperature at which weathering becomes active will de- 
termine the surface temperature by adjusting the atmospheric greenhouse CO2 
level for the long-wave outgoing radiation balancing the incoming solar radi- 
ation. The thermostat will regulate temperature until the solar flux is so 
large that it can maintain a surface temperature above the temperature where 
weathering becomes active without the CO2 greenhouse warming. From that 
point in time the surface temperature will increase with the increasing solar 
flux. The greenhouse thermostat is keeping the climate warmer than the sta- 
ble Snowball Earth situation where the ice albedo keeps the surface globally 
below the freezing point of water. Relatively constant temperatures over ge- 



ological time determined by geochemical thermostats could be important for 
initiation of biological life and potentially widen the habitable zone around a 
sun-like star. When bacterial type life or the like is established on a planet 
it could provide a thermostat by itself, as suggested by the gaia hypothesis 
(Lovelock fc Margulis, 1982). Life on Earth is today a dominant player in 



controlling the atmospheric levels of O2, CO2 and CH4. 
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FIGURE CAPTIONS 

Fig. 1 The incoming (red) and outgoing (blue and green) fluxes as a function 
of global temperature. When the greenhouse warming is below 22 K 
the warm stable climate state at T c disappears through a saddle-node 
bifurcation and the 'deep freeze' state at T a is the only stable climate 
state. 

Fig. 2 The bifurcation diagram for the energy balance. If the greenhouse warm- 
ing fall below Tg the climate will fall into the Snowball Earth climate. 

(2) 

The greenhouse warming has to exceed Tg for the planet to leave the 
Snowball Earth climate. The dodded arrows indicate a hysteresis loop. 

Fig. 3 The stability of the warm climate state against random fluctuations is de- 
termined by the radiative cooling feedback represented by the parameter 
—a. 

Fig. 4 The surface temperature has, except for a few Snowball Earth episodes in 
the late Neoproterozoic, been above the freezing point of water for most 
of Earths geological history since the heavy bombardment epoch. The 
surface temperature has been governed by the greenhouse thermostat. 
The thermostat was operational until the solar luminosity became strong 
enough perhaps 1.5-1 Ga BP. The present atmospheric content of CO2 
and other gasses is regulated by biology itself, not included in the graph. 
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